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Three consecutive plaque purifications of four chimeric yellow fever virus-dengue virus (ChimeriVax-DEN) 
vaccine candidates against dengue virus types 1 to 4 were performed. The genome of each candidate was 
sequenced by the consensus approach after plaque purification and additional passages in cell culture. Our 
data suggest that the nucleotide sequence error rate for SP6 RNA polymerase used in the in vitro transcription 
step to initiate virus replication was as high as 1.34 x 10~ 4 per copied nucleotide and that the error rate of 
the yellow fever virus RNA polymerase employed by the chimeras for genome replication in infected cells was 
as low as 1.9 x 10 ~ 7 to 23 x 10 ~ 7 . Clustering of beneficial mutations that accumulated after multiple virus 
passages suggests that the N -terminal part of the prM protein, a specific site in the middle of the E protein, 
and the NS4B protein may be essential for nucleocapsid-envelope interaction during flavivirus assembly. 



A fundamental feature of RNA viruses is the error-prone 
nature of replication of their genomes, resulting in adaptation 
and rapid evolution. Various in vitro, in vivo, and ex vivo 
methods have been used to estimate mutation rates of viral 
RNA-dependent RNA polymerases. The generally accepted 
average value of 10 ~ 4 to 10~ 5 mutations per nucleotide (nt) 
per round of RNA replication (reviewed in references 6, 7, 23, 
and 28) results in substitution of 0.1 to 1 nt per each synthe- 
sized genomic RNA molecule of —10 kb. Mutation rates vary 
between 10~ 3 and 10~ 6 in individual studies, depending on the 
virus and method used, although none of these methods pro- 
vide an exact measure of the error rate. In vitro methods that 
utilize purified viral RNA polymerases may overestimate the 
error rate due to suboptimal in vitro enzymatic conditions and 
experimental complexity. In ex vivo methods, such as sequenc- 
ing of phenotypic mutants generated in cell culture and anal- 
ysis of monoclonal antibody escape mutants or of rates of 
reversion of phenotypic markers, the number of rounds of 
RNA replication in infected cells is not precisely known. Meth- 
ods that involve a reverse transcription-PCR step with cloning 
of cDNA in bacteria prior to sequencing can be misleading 
because of the accumulation of additional mutations intro- 
duced by reverse transcription-PCR. 

The Flavivirus genus of the Flaviviridae family consists of 
about 70 viruses of which 40 have been associated with human 
illness (reviewed in reference 2). Flaviviruses are small envel- 
oped plus-strand RNA viruses (reviewed in reference 18). The 
viral particle contains a nucleocapsid composed of viral RNA 
and capsid protein C, surrounded by a lipid bilayer in which the 
envelope protein E and membrane protein M are embedded. 
The genomic RNA is about 11,000 nt long and contains a 
single long open reading frame (ORF). Translation of the 
ORF produces a single polyprotein precursor containing viral 
proteins in the order: C-prM-E-NSl-NS2A-NS2B-NS3-NS4A- 
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2K-NS4B-NS5, where C through E are the structural compo- 
nents of the virion (prM is a precursor for M). Nonstructural 
(NS) proteins NS1 through NS5 are required for virus repli- 
cation in the cytoplasm of infected cells (18). Immature virions 
are produced by budding into the lumen of the endoplasmic 
reticulum (ER) and movement of virus to the cell surface 
through the exocytosis pathway. The structure of the flavivirus 
particle has been resolved by cryoelectron microscopy and 
fitting the known structure of the E protein (29) into the 
electron density map (15). 

The ChimeriVax technology has been used to create atten- 
uated flavivirus vaccines using the yellow fever virus (YF) 17D 
vaccine virus as a vector in which the prM-E envelope protein 
genes are replaced with those from a heterologous flavivirus (4, 
10, 11, 24-26). To produce a chimeric virus, the viral RNA 
genome is first synthesized by in vitro transcription of an ap- 
propriately engineered DNA template with bacteriophage SP6 
RNA polymerase. Following transfection of cells with the in 
vitro RNA transcripts, virus replication begins with genome 
amplification by the YF 17D RNA-dependent RNA polymer- 
ase (YFpol). In this study, yellow fever virus-dengue virus 
(ChimeriVax-DEN) type 1 to 4 chimeras were prepared by 
transfecting Vero cells with in vitro RNA transcripts. The 
resulting viruses were cloned by three consecutive plaque pu- 
rifications. Both cloned and uncloned viruses were sequentially 
passaged to evaluate their genetic stability during prolonged 
growth in cell culture. Each virus was sequenced by the con- 
sensus method several times at different passages, including 
before and after plaque purification of the cloned variants. We 
observed that YFpol introduced very few mistakes in the ge- 
nome during plaque purification steps, indicating a high fidelity 
of the polymerase. Characterization of genomes with YFpol 
mistakes and mutations that accumulated during multiple pas- 
sages shed new light on the mechanism of flavivirus evolution, 
RNA replication, and particle assembly. Data from these stud- 
ies challenge the widely held concept that RNA genomes ac- 
cumulate substitution of up to 1 nt for each synthesized 
genomic RNA molecule. 
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FIG. 1. Passages of the cloned and uncloned Chimed Vax- 
DEN^ooo viruses. PI virus was obtained by electroporation of Vero 
cells with in vitro RNA transcripts synthesized using SP6 RNA poly- 
merase. Uncloned P2 virus was generated by one additional passage. 
Cloned PMS candidates (clones A and B) were produced by three 
consecutive plaque purifications starting from PI virus followed by two 
additional passages. Genetic stabilities of these cloned and uncloned 
variants were examined in additional downstream virus passages. Chi- 
meriVax-DEN22ooi, -DEN3, -DEN4, and -DEN1 viruses were pas- 
saged similarly, except that plaque purification passages (P3 to P5) 
were initiated from the uncloned P2 viruses. 



Virus passages* A diagram of passages used to prepare un- 
cloned passage 2 (P2) and two cloned P6 Pre-Master Seed 
(PMS) variants of ChimeriVax-DEN2 2000 virus is shown in 
Fig. 1. Uncloned and cloned ChimeriVax-DEN2 2001 , -DEN3, 
-DEN4, and -DEN1 PMS viruses were prepared similarly, ex- 
cept that plaque purification passages (P3 to P5) were initiated 
at the uncloned P2 level rather than PI. Vero cells provided by 
Aventis Pasteur (Lyon, France) were used in all experiments. 
Cells and viruses were grown in Earl's minimum essential me- 
dium (EMEM) supplemented with L-glutamine, MEM nones- 
sential amino acids (Gibco), and 10% fetal bovine serum (FBS; 
HyClone). Medium 199 (Gibco) containing 3% FBS was used 
for plaque purifications done in six-well culture plates by the 
agarose double overlay method with neutral red staining de- 
scribed previously (24). 

In vitro RNA transcripts were synthesized with SP6 RNA 
polymerase (Epicentre) on full-length DNA templates for the 
four ChimeriVax-DEN viruses prepared using two- or three- 
fragment ligations (11). No PCR steps were involved in prep- 
aration of the templates, and each plasmid was grown from a 
single bacterial colony. The PI viruses were obtained by elec- 
troporation of 1.5 X 10 7 cells (which had been trypsinized and 
washed with phosphate-buffered saline) with 50 to 200 ng of 
the in vitro RNA transcript, followed by seeding of the cells 
into T75 or T25 flasks. PI virus was harvested on days 3 to 5 
posttransfection. The uncloned P2 PMS ChimeriVax-DENl to 
-DEN4 viruses, cloned P6 ChimeriVax-DEN2 2000 virus, and 
cloned P7 ChimeriVax-DENl, -DEN3, -DEN4, and -DEN^^ 
PMS viruses were obtained by infection of cells at a controlled 
multiplicity of infection (MOI) of 0.001 PFU/cell. Titers of 



these virus samples were in excess of 10 6 PFU/ml. Up to 13 
additional (genetic stability) passages of these viruses were 
done without intermediate virus titrations at an assumed MOI 
of 0.001 PFU/cell. All virus harvests were collected when ~5 to 
10% cytopathic effect was observed. Virus-containing media 
were clarified by ]ow-speed centrifugation, diluted with an 
equal volume of FBS, and stored at -70°C. Up to 10 clones 
were prepared for each chimera by three sequential plaque 
purifications. Plaque picks (agarose plugs) were placed in 0.4 
ml of medium 199 (Gibco) containing 50% FBS, vortexed for 
several seconds, and frozen at -70°C. To sequence viruses 
obtained after the first plaque purification, 100 uJ of the indi- 
cated plaque pick supernatant (virus eluted from agarose plug) 
was used to infect cells in T25 flasks, and the amplified viruses 
were harvested and sequenced. Otherwise, all samples of virus 
passages shown in Table 1 were sequenced without additional 
expansion. 

Sequence analysis. Consensus sequencing of virus samples 
was performed as previously described (27a). Sequence data 
were aligned and analyzed using Sequencher 3.0 (GeneCodes) 
software. Nucleotide heterogeneities were registered when a 
heterogeneous signal was observed in all chromatograms rep- 
resenting both plus and minus cDNA strands. For some viruses 
(Table 1), only the first of the five cDNA amplicons (fragment 
I) was prepared and sequenced, including the 5' untranslated 
region (5' UTR) and the structural protein genes. In several 
cases, such as viruses obtained after the first plaque purifica- 
tion or some intermediate genetic stability passages, only the 
regions where mutations were detected at later passages were 
sequenced, assuming that no other mutations were present in 
the rest of the genome (or fragment I). 

Consensus full-genome sequencing, which does not detect 
random mutations, of the uncloned P2 viruses revealed no 
unexpected mutations (Table 1). Mutations accumulated in the 
uncloned ChimeriVax-DENl to -DEN3 variants during their 
further passages to study genetic stability (see passages after 
P2 in Table 1). 

Since no mutations were detected in the uncloned P2 Chi- 
meriVax-DEN2 2000 virus, no mutations should be detected in 
the consensus sequence of PI virus used to start plaque puri- 
fications. Nevertheless, each of the cloned A and B P6 Chi- 
meriVax-DEN2 20O0 viruses had a 1-nt change: G to T at nt 
7734, resulting in an R-to-M amino acid change at NS5-29 
located in a relatively conserved N-terminal region of NS5 and 
a silent A-to-T change at nt 6643, respectively (Table 1). Two 
additional genetic stability passages of clone A virus resulted in 
accumulation of subpopulations containing changes at nt 515 
and 2135 detected at P8. Changes in clone A had potentially 
significant phenotypic implications. Therefore, further work 
with this clone was terminated. In clone B virus at P9, an 
additional G/C heterogeneity was observed at nt 517. The 
G-to-C nt change at this position resulted in an M-to-I amino 
acid change at prM-12 in the N-terminal part of prM that is 
cleaved during virus maturation. The effect of this mutation is 
similar to that of the nt 515 change in clone A, resulting in a 
similar mutation (M to L) at the same residue of prM. An 
M-to-L change at this residue was also seen in uncloned P15 
ChimeriVax-DEN2 2O01 virus (Table 1). This heterogeneity be- 
came a dominating change in virus at P14. An L-to-F change at 
amino acid 177 of NS4B protein was selected during virus 
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TABLE 1. Nucleotide and amino acid changes in uncloned and cloned ChimeriVax-DEN viruses 



Virus 



Version 



Mutations detected at sequenced passages 0 



Uncloned P2, none*; P5, none*; P10, A 17J6 T (K^^M), G, 287 A (V NSI . 275 M)* 

Clone A P2, G 7734 T (R NS5 . 29 M) C ; P6, G 7734 T (R NS5 . 39 M) b ; P8, A 515 /C (M prM _ 12 /L), CWG (QbWE), 

G 7734 T (R N S5-29M) & 

Clone B P2, A^T (silent)*; P6, A 6643 T (silent) 6 ; P9, G 517 /C (M prM . 12 /I), A^T (silent)*; P14, G 517 C 

(M prM -i 2 I)> A^T ^silent), C, 477 T (L NS4B . 177 F) C ; P18, G 517 C (M prM . 12 I), W (silent), 

Uncloned P2, none*; P15, A 5l5 C (M M . 12 L), G 132 jA (silent), A 1590 G (K E _ 204 R), A 4046 T (I NS2 a-i76 L )» 

T 5485 C (silent), C, 888 T (A 2K _ 20 V), C 7427 T (L^^F), G^/T (silent), C 103O 7T (silent) 6 

Clone A P6, none 6 ; P7, none*; P10, T 5S1 G (1^-24 V), G 1730 T (V^F) 6 ; P20, T 551 G (L prM _ 24 V), 

G a730 T (V^F)* 

Clone B P3, none c ; P6, A 1030 G (silent), C 1789 T (silent) rf ; P15, A 1030 G (silent), A, 590 G (K^^R), 

(silent), T 2444 C (silent), C 7443 T (A NS4B . 182 V), A 8677 G (silent) 6 
Uncloned P2, none 6 ; P5, none c ; P10, C 564 A (A-m^D), T^C (Y NS1 _ 35 H) C ; P15, C 564 A (A^^D), 

(Y NS1 - 35 H), A« 502 /G (silent), A 7231 /G (Wb-u-ZM)* 
Clone A P3, C^T (silent)'; P6, C^T (silent) 6 ; P7, C^T (silent) 6 ; P10, C^T (silent) c ; P15, 

C^T (silent), (WT (silent)**; P20, Q^T (silent), C 7421 T (Ln$4b-i77F)> CWT 
(silent) 6 

Uncloned P2, none 6 ; P15, none 6 

Clone B P6, none 6 ; P7, none 6 ; P10, none 6 ; P20, none 6 

Uncloned P2, none 6 ; P5, A 1590 /G (K E . 204 /R), G 1730 /T (Ve.^/F)*; PI 5, A 159 qG 

G l73( /r (v E . 

251/F), A 7237 /G (I N S4B-113/M), 6,4^ (P NS 4B-WS) 6 

Clone A P3, G 1730 T (V^F), C^A (L^I)'; P7, G 1730 T (V E . 25l F), C 2282 A (L^I)* 

Clone B, D P3, G 1730 T (Ve. 251 F) c ; P7 or P6, G 1730 T (V &2S1 F)* 

Clone C P3, G 1912 T (E^D), G 2O30 T (V M51 L) C ; P6, G 19l2 T (E E . 3U D), G 203O T (V^ S1 L) 6 

Clone E P3, A 1590 G (K^R), A 3952 T (siient) c ; P6, A 1590 G (K MM R), A 3952 T (silent) 6 

Clone F P3, C^T (silent), A, S90 G (K E . 204 R) d 

Clone G P3, G 1730 T (V^F)* 

Clone H P3, G )912 T (E^D), G 2030 T <V M51 L) d 

Clone I P3, A 1590 G (K E . 204 R) rf 

Clone J P3, none 6 ; P6, none 6 ; P7, none 6 ; P10, A 1590 G (K E . 204 R) 6 ; P20, A 1590 G (K E . 204 R), G 6966 /T 
(Sns4b-23/I)> G 7190 /A (V NS4B . 98 /I) 6 



ChimeriVax-DEN2 2O0O 
ChimeriVax-DEN2 2001 

ChimeriVax-DEN3 

GhimeriVax-DEN4 
ChimeriVax-DENl 



; passage number is shown first; any rautation(s) detected follow the passage number. Nucleotide changes at indicated positions (numbering from the beginning 
genome) and amino acid changes (in parentheses) at indicated protein residues [from the N terminus of each protein; amino acid numbering of YF proteins 



"The 1 
of the genome) i 

as in references 3 and 17] are shown. A heterogeneous signal is indicated by a slash. Data for clones after the first and/or last plaque puriBcations are in bold type. 
b Full genome sequenced. 

c Sites of mutations found after last plaque purification (in samples after first plaque purification) or at a later genetic stability passage were sequenced only. It can 
be assumed that the full genome (or fragment I) sequence is known. 
d Only fragment I sequenced, including the entire structural protein coding region. 

c In independent passages of the DEN3 clone A virus to P15 from a Master Seed (large-scale manufacturing P8), accumulation of A 1584 G (K^^R; identical to the 
E-204 mutation in the DEN1 and DEN2 chimeras) and C7 437 T (A NS4B -i82 v ) mutations was documented, in addition to the silent nt 6607 change (F. Mitchell, data 
not shown). 



passage to P14. The P18 of clone B had a consensus sequence 
identical to that of P14. It was significant that the P2 samples 
of both clones A and B harvested after the first plaque purifi- 
cation contained the same mutations as the P6 viruses at nt 
7734 and 6643, respectively (Table 1). 

Starting plaque purification with the uncloned P2 virus, 
rather than with PI virus, immediately yielded acceptable 
cloned vaccine candidates of ChimeriVax-DEN2 2001 (clone A), 
-DEN3 (clone A), and -DEN4 (clone B; clone A was not 
sequenced) (Table 1). Only clone A of the ChimeriVax-DEN3 
chimera had one silent nucleotide change at position 6607 
detected at both P3 after the first plaque purification and P6 
after the third plaque purification. In contrast, 10 clones of 
ChimeriVax-DENl were examined to find one without unac- 
ceptable mutations. Only clone J of the ChimeriVax-DENl 
chimera had no mutations at both P3 and P6. The entire 
genomes of clones A, B, C, and E were sequenced at P6, and 
the presence of mutations was checked at P3. Only the struc- 
tural protein coding region of clone D was sequenced at both 
P3 and P6. In all cases, mutations present at P6 were present 



in P3 viruses (Table 1). (Fragment I of clones F to I was 
sequenced only at P3.) 

Some clones of the ChimeriVax-DENl chimera had identi- 
cal mutations, such as the E-251 V-to-F change in clones A, B, 
D, and G; the E-204 K-to-R change in clones E, F, and I; and 
the E-311 E-to-D and E-351 V-to-L changes in clones C and H. 
These changes were not detected in the consensus sequence of 
the starting uncloned P2 virus. Mutations E-204 and E-251 
were first detected as heterogeneities in the uncloned P5 virus. 
The E-204 mutation became dominant in the uncloned P15 
virus, and it also appeared in clone J virus by P10 (Table 1). 
The rate of accumulation of the E-204 mutation was studied in 
more detail in additional passaging experiments using clone J. 
These studies indicated this mutation was infrequently intro- 
duced by YFpol (C. Penal and F, Mitchell, data not shown). 
Since the E-204 mutation and all other mutations found in 
clones A to I at P3 appeared as clear mutations, they must have 
been present in individual particles of the uncloned P2 virus 
population before the first plaque purification. The fact that 
the E-204, E-251, E-311, and E-351 mutations were not de- 



Vol. .78, 2004 



NOTES 1035 



tected in the consensus P2 sequence despite the high propor- 
tion of these mutations in the plaques picked was most likely 
due to the low sensitivity of consensus sequencing. Our prior 
experience with quantitation of mutant subpopulations using 
the MAPREC method (5) demonstrated that the level of de- 
tection of mutations by sequencing is approximately 10 to 15%. 
During the first plaque purification step, a few slightly larger 
plaques were picked together with smaller plaques which could 
account for some of the observed variation. The possibility that 
any of these mutations originated in plasmid templates is un- 
likely because of the high fidelity of bacterial DNA synthesis 
(7, 23). Ten individual plasmid clones from the ChimeriVax- 
DEN1 DNA template were prepared in Escherichia coli and 
sequenced. None of the clones contained the mutations de- 
tected in clones A to I of the virus or any other mutations (M. 
Parsons, data not shown). 

Clone B of ChimeriVax-DEN22 001 was prepared and se- 
quenced at P3 after the first plaque purification, P6 after the 
third plaque purification, and P15. The P6 virus had two silent 
nucleotide changes at nt 1030 and 1789 which were not present 
at P3 (Table 1). These mutations occurred during plaque pu- 
rification passages, and thus represented mistakes of YFpoI. 
These were the only two mistakes introduced by YFpol during 
three plaque purifications of all the cloned chimeras described. 

( ' The cloned and uncloned ChimeriVax-DEN4 viruses did not 
accumulate mutations during genetic stability passages. In con- 
trast, both cloned and uncloned variants of the other chimeras 
accumulated various mutations in genetic stability passages. 
Most of these mutations resulted in amino acid changes, indi- 
cating that virus adaptation is driven by forces operating at the 
protein function level. Some of the mutations observed among 
different chimeras coincided (Table 1). For instance, the same 
E-204 mutation accumulated in uncloned and clone B Chi- 
meriVax-DEN2 2O01 (detected at P15 of both variants), un- 
cloned and clone J of ChimeriVax-DENl (at P5 and P10, 
respectively), and ChimeriVax-DEN3 (the E-202 mutation; see 
footnote e to Table 1). This suggests convergent mechanisms of 
adaptation of the ChimeriVax-DENl to -DEN3 chimeras to 
Vero cells. Many diverse mutations accumulated in uncloned 
and cloned versions of the same chimera or in different clones 
of a chimera (Table 1). The latter observation indicates that 
one virus can evolve in different directions under similar 
growth conditions. It is possible that several different muta- 
tions can improve virus replication and that further virus evo- 
lution depends on which mutation occurs first. 

: Classification of mutations. The interpretation of mutations 
observed at different passages of cloned and uncloned Chi- 
meriVax-DEN viruses is summarized in Table 2. No mutations 
were detected in the consensus sequences of all chimeras in 
uhcloned P2 virus. As a rule, mutations in cloned viruses de- 
tected after three consecutive plaque purifications (at P6 or 
P7), were also observed after the first plaque purification. We 
believe that these mutations were present in individual virions 
present in the uncloned PI or P2 virus population obtained 
after electroporation of cells with in vitro RNA transcripts. 
These mutant viruses were randomly selected by the first 
plaque purification, as no additional mutations occurred dur- 
ing subsequent plaque purification passages. Since each act of 
plaque formation requires multiple amplifications of the 
genomic RNA in the initiating viral particle, YFpoI must have 



TABLE 2. Assigned classification of mutations detected during 
preparation of ChimeriVax-DEN viruses and the genetic stability 
passages shown in Table 1 



Mutation detected 



Type of mutation 



In uncloned 
P2 virus 



After the 
first plaque 
purification 0 



After the In genetic 
third plaque stability 
purification 6 passages 



Introduced by SP6 

polymerase (17 

nucleotide 

changes in 11.5 

genomes) 
Introduced by 

YFpol (2 

nucleotide 

changes in 10.5 

genomes) 
Beneficial mutation* 

(15 and 19 

nucleotide 

changes in the 

cloned and 

uncloned variants, 

respectively) 



" P2 plaque purification step of cloned ChimeriVax-DEN22ooo or of cloned 
ChimeriVax-DENl, -DEN2 2001 , -DEN3, and -DEN4 (Table 1) (Fig. 1). 

b P6 of cloned ChimeriVax-DEN2 20 oo or P6 or P7 of cloned variants of other 
chimeras (Table 1). 

c Beneficial mutations that accumulated in genetic stability passages of cloned 
variants were introduced by YFpol, while it is unknown which polymerase (SP6 
or YFpol) introduced mutations that accumulated in uncloned variants. The 
E-204 and E-251 mutations beneficial for ChimeriVax-DENl were also found 
after plaque purifications in some of the clones of the chimera. 



high fidelity. Each act of plaque formation is a very harsh virus 
passage in which replication is allowed to proceed until com- 
plete destruction of participating cells. Therefore, it is unlikely 
that mutations were introduced by YFpol in the process of 
virus replication prior to the first plaque purification. Most 
likely, they were introduced by SP6 RNA polymerase during in 
vitro RNA synthesis. This rule applies to all cloned viruses, 
except for clone B of ChimeriVax-DEN2 2001 (Table 1). This 
virus had two silent nucleotide changes at P6 level (nt 1030 and 
1789), which were not present at P3. Thus, these two mutations 
were the only detected mistakes of YFpol that occurred during 
plaque purifications. 

All mutations in individual cloned viruses selected by plaque 
purifications were stably fixed. As a rule, these mutations did 
not accumulate in other clones or uncloned variants during 
genetic stability passages. Other mutations accumulated in ge- 
netic stability passages of most clones and uncloned variants. 
Subpopulations containing such mutations outgrew the non- 
mutant viruses in the population, and therefore such mutations 
are beneficial. All the chimeras are viable immediately after 
transfection. However, analyses of some of the accumulating 
virus variants showed that these variants produce larger 
plaques in Vero cells and higher peak titers compared to non- 
mutants (by up to 1 log 10 PFU/ml; data not shown). Thus, 
there is room for improvement of growth of these artificial 
chimeras. 

Some mutations in ChimeriVax-DENl clones A to I de- 
tected after the first and the third plaque purification (Table 1) 
did not clearly fit the classification in Table 2. The E-204 and 
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E-251 mutations were detected in more than one clone and 
also accumulated in genetic stability passages of the uncloned 
virus. The E-204 mutation also accumulated during continuous 
passage of ChimeriVax-DENl clone J. Thus, these mutations 
are beneficial for the ChimeriVax-DENl chimera replication. 
Since mutations are introduced by YFpol infrequently, these 
changes were likely mistakes of SP6 polymerase amplified 
prior to the first plaque purification. 

Of the 17 nucleotide changes considered to be mistakes of 
SP6 polymerase, 5 were transitions and 12 were transversions. 
Thirteen of the mutations were changes from A, G, or C 
nucleotides to T, with a G-to-T change being the most frequent 
(nine mutations). This suggests that a misincorporation of T 
was the most common mistake during in vitro transcription. Of 
the 17 nucleotide changes introduced by YFpol, including the 
two mistakes that occurred during plaque purifications plus 
beneficial mutations that accumulated during genetic stability 
passages of the clones, 11 changes were transitions and 6 were 
transversions. The most frequent changes were A to G (four 
mutations) and C to T (five mutations), which theoretically 
should be the two most common mistakes during late plus- 
strand RNA synthesis due to T • G and G * T mispairing, 
respectively. This observation may point toward the existence 
of distinct mechanisms for synthesis of double-stranded RNA 
(with high fidelity), representing most rounds of amplification 
in infected cells, and late transcription (with low fidelity). The 
majority of rounds of viral RNA synthesis should occur during 
the latent phase of infection accompanied by exponential du- 
plication of the double-stranded replicative-form (RF) RNA. 
During the late linear phase, RNA synthesis is mostly single- 
round transcription on existing minus-strand template that is 
stable for at least 20 h and is renewed slowly (reviewed in 
references la and 30). 

Estimates of error rates for SP6 RNA polymerase and 
YFpol. Sequences of 11.5 equivalents of complete virus ge- 
nome determined after the first plaque purification (sequences 
of fragment I are counted as 0.25 genome) contain a total of 17 
nucleotide changes (Table 1) considered to be mistakes of SP6 
polymerase introduced during a single round of DNA-depen- 
dent RNA synthesis in vitro. Thus, the error rate of SP6 RNA 
polymerase is 17/11.5 = 1.48 misincorporations per one 
genomic RNA molecule synthesized (—11,000 nt long), or 1.34 
x 10~ 4 per nucleotide copied, which is comparable to the 
reported error rate of T7 RNA polymerase (0.5 X 10" 4 [13]). 
However, it could be somewhat lower because all of the mu- 
tations may not be due to mistakes of SP6 polymerase and the 
accumulation of mutations in ChimeriVax-DENl may have 
occurred prior to the first plaque purification. On the other 
hand, if plaque purifications of ChimeriVax-DEN2 2001 , 
-E>EN3, and -DEN4 had been done at PI, rather than P2, more 
mistakes of SP6 polymerase may have been detected. 

Two nucleotide changes were present in a total of 10.5 full 
genome equivalents that occurred during plaque purifications 
and therefore were clearly introduced by YFpol (Table 1). To 
estimate the YFpol error rate, we need to know how many 
times a genomic RNA molecule is copied during one act of 
plaque formation. Chimeras grow to titers of ~10 7 PFU/ml in 
Vero cells; therefore, a plaque pick of roughly 100 jxl in volume 
will contain ~10 6 infectious particles. We assume that 100 to 
1,000 times more plus strands are synthesized (10 8 to 10 9 mol- 

r 



J. Virol. 

ecules) and are not accounted for in infectious virus because of 
the loss of infectivity of mature secreted particles, as well as the 
accumulation of immature viral particles and possible RNA 
turnover inside the cells (reviewed in references la and 30). If 
most rounds of RNA amplification are semiconservative du- 
plications of the RF RNA, as many as 25 to 30 rounds of 
synthesis can take place (e.g., 10 9 molecules ** 2 30 , thus re- 
quiring 30 duplications) during one plaque formation, in ~10 3 
to 10 4 participating cells. In support of our estimate of the total 
number of plus strands synthesized, 20 to 40 jxg of protein E 
per 10 7 cells is produced following infection with tick-borne 
encephalitis virus (1, 14a). Since each viral particle contains 
180 protein E monomers (15), up to 2 X 10 5 viral particles 
(equal to secreted plus strands) are produced by each cell. This 
number does not include all plus strands synthesized in in- 
fected cell. 

By dividing the number of YFpol mistakes (two mistakes) by 
the number of sequenced full genome equivalents (10.5 ge- 
nomes), the estimated number of rounds of RNA synthesis 
during one plaque formation (25 to 30 rounds), and the num- 
ber of plaque purification steps for each clone (three steps), 
the resulting error rate of YFpol is 2/[10.5 x (25 to 30) x 3] = 
0.0021 - 0.0025 misincorporation per genomic RNA molecule 
(11,000 nt long) per synthesis, or 1.9 x 10~ 7 to 2.3 X 10" 7 per 
copied nucleotide. With our system, we can detect only non- 
lethal mutations. We do not believe, however, that lethal mu- 
tations could constitute a percentage significant enough to 
affect the estimate, considering the overall low occurrence of 
mutations during plaque purifications. It is known that the 
distribution of mutations in RNA viruses is shifted toward 
neutral or slightly deleterious mutations rather than strongly 
deleterious or lethal mutations that are quickly eliminated by 
competition during virus replication; besides, a small fraction 
of mutations is beneficial (8, 9). To conclude that there could 
be a large proportion of undetected lethal mutations, we must 
observe many synonymous mutations, such as those that do not 
change the amino acid encoded, which is not the case. 

Our data indicate that the mutation rate of YFpol is on the 
low end of the spectrum among RNA polymerases. This ob- 
servation is in accord with previous studies demonstrating the 
high genetic stability of YF (12, 27, 33) and makes YF 17D an 
excellent vector for the construction of recombinant vaccines. 
It is believed that the evolution of YF and other flaviviruses in 
nature is constrained by the need of alternation of vertebrate 
and arthropod hosts. Our data extend this hypothesis further, 
suggesting that in the course of flavivirus evolution, variants 
with high fidelity of their RNA polymerases were selected to 
facilitate the necessary continuous host alternation. 

Hot spots for beneficial amino acid changes. All beneficial 
amino acid changes that accumulated during genetic stability 
passages of the cloned and uncloned viruses were plotted in 
Fig. 2. A conspicuous clustering of mutations was observed, 
specifically in the N-terminal part of prM in ChimeriVax- 
DEN2 and -DEN3 chimeras, the central part of the E protein 
in ChimeriVax-DENl, -DEN2, and -DEN3, and in the NS4B 
protein in ChimeriVax-DENl, -DEN2, and -DEN3. One could 
argue that if we adapted wild-type DEN viruses or YF 17D to 
Vero cells, the same mutations could be expected. This ap- 
pears not to be the case. Specific mutations in the YF-specific 
NS4B protein did not accumulate in the ChimeriVax-DEN4 
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FIG. 2. Clustering of beneficial amino acid changes in the four ChimeriVax-DEN chimeras that accumulated in all cloned and uncioned viruses 
during genetic stability passages. The locations of mutations in the polyprotein and genomic RNA in each of the chimeras were plotted according 
to their positions in the genome. 



virus. A recent study on adaptation of tick-borne encephalitis 
virus grown in BHK-21 cells demonstrated that the majority of 
adaptive mutations resulted in an increase in the net positive 
charge of the E protein, presumably facilitating the protein 
interaction with glycosaminoglycans on the cell surface (22). 
None of the observed mutations in the E protein in Chi- 
meriVax-DENl (E-204 K-to-R and E-251 V-to-F mutations), 
ChimeriVax-DEN2 (E-204 K-to-R, E-246 K-to-M, E-251 V- 
to-F, and E-386 Q-to-E mutations), and ChimeriVax-DEN3 
(E-202 K-to-R mutation) were of this kind. Some of the mu- 
tations could affect the ability of the viruses to antagonize 
cell-specific (e.g., RNA interference) or -nonspecific antiviral 
mechanisms. It is still unclear why no mutations accumulated 
in the ChimeriVax-DEN4 chimera. 

We favor another possibility — that these mutations are ad- 
aptations in response to replacing the YF-specific envelope 
with the heterologous DEN envelopes. An imperfect match 
between the YF-specific capsid and the DEN1- to DEN3- 
specific envelopes dictated adaptation during virus growth. 
The DEN4 chimera needed no adaptation due to an incidental 
perfect match between the capsid and DEN4 envelope. The- 
oretically, this is possible, despite the fact that DEN4 is a type 
of the same virus species because the specific requirements for 
interaction between the capsid and envelope are unknown, 
j If the latter hypothesis is correct, it is interesting that no 
mutations accumulated in the YF-specific capsid protein C. 
The mutations in the N-terminal part of prM could have func- 
tional significance in light of recent data indicating that prM 
functions as a chaperone necessary for proper folding of E 
(20). These mutations could facilitate a more appropriate fold- 
ing of E such that the envelope more precisely conforms to the 
architecture of the nucleocapsid. All mutations in the cluster 
located in the middle of E protein are on the outward surface 
of the protein in the dimerization domain II (29). The effects 
of these changes on the structure and function of the E pro- 
teins of chimeric viruses need further investigation. 

. Most intriguing is the cluster of mutations in NS4B (Fig. 2), 
because it suggests that this protein plays a role in virus as- 
sembly. Flavivirus replication induces proliferation of smooth 
ER membranes and the accumulation of characteristic mem- 



branous structures described as convoluted membranes (CM), 
paracrystalline arrays (PC), and vesicle packets (VP) (21, 31, 
32). Viral double-stranded RNA together with NS1, NS3, NS5, 
NS2A, and NS4A proteins colocalized with VP structures, in- 
dicating that they are the main sites of RNA synthesis, while 
CM and PC structures were proposed to be main sites of 
proteolytic cleavage (14, 21, 32). NS4B appears to have a 
distinct distribution, associated with ER membranes through- 
out the cytoplasm (31). Recently, NS2A was implicated in virus 
assembly, because specific mutations in this protein abolished 
production of infectious particles (16). In addition, it was pro- 
posed that NS3 also plays a role in assembly (14, 16, 19). Since 
NS2A and NS3 are both located in the sites of RNA synthesis 
and taking into account that both these proteins bind viral 
genomic RNA (21), it is tempting to speculate that NS2A and 
NS3 are responsible for genome encapsidation, while NS4B 
comes into play during the next step, budding, possibly by facili- 
tating specific interaction of the nucleocapsid and envelope, 
NS4B could also facilitate particle formation indirectly, via its 
involvement in the formation of virus-induced membrane struc- 
tures essential for the secretion pathways (14b, 20a). 

Finally, the suggestive nature of these sequencing data 
should be recognized. The evolutionary, enzymological, and 
molecular aspects discussed above require further investiga- 
tion. The process of accumulation of viral RNAs in infected 
cells needs to be precisely quantitated using modern ap- 
proaches such as real-time PCR. This could establish the exact 
number of rounds of RNA synthesis, confirm or correct our 
estimate of YFpol fidelity, and add to a better understanding 
of the molecular mechanisms of flavivirus RNA synthesis. The 
high fidelity of YFpol requires confirmation by other available 
methods, including using wild-type YF. The hypothetical role 
of NS4B (as well as the elements in prM-E) in particle assem- 
bly requires direct demonstration. 

This work was funded by Aventis Pasteur, Lyon, France. 
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